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The critical phenomena in strongly interaction matter are generally in- 
vestigated using the mean-field model and are characterized by well defined 
critical exponents. However, such models provide only average properties 
of the corresponding order parameters and neglect altogether their possible 
fluctuations. Also the possible long range effect are neglected in the mean 
field approach. Here we investigate the critical behavior in the nonexten- 
sive version of the Nambu Jona-Lasinio model (NJL). It allows to account 
for such effects in a phenomenological way by means of a single parameter 
q, the nonextensivity parameter. In particular, we show how the nonexten- 
sive statistics influence the region of the critical temperature and chemical 
potential in the NJL mean field approach. 

PACS numbers: 21.65.+f; 26.60.+c; 25.75.-q 

1. Introduction 

Critical phenomena in strongly interaction matter are of great interest 
nowadays, cf., for example, [HI2]- They are usually described by a mean field 
type of theories^. Such theories are based on the usual Boltzmann-Gibbs 
(BG) statistical mechanics and reflect only behavior of the mean field, i.e., 
are not able to accommodate effects of possible fluctuations and/or cor- 
relations caused, among others, by the smallness of the sample of matter 
under consideration, by its rapid evolution, or by limitations of the avail- 
able phase space. All these factors (both separately and taken together) 
render the spatial configuration of the system being far from uniform and 
prevent the global equilibrium from being established. Nevertheless, it is 
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know in the literature that it is possible to maintain simplicity of the statis- 
tical description and, at he same time, to account for these effects provided 
one uses a nonextensive version of the statistical mechanics, for example 
the one proposed by Tsallis [5j (which we shall use here). In this ap- 
proach action of all the above mentioned factors is summarily accounted 
for by one additional parameter q, the nonextensivity parameter. It does 
not differentiate between the particular dynamical phenomena responsible 
for departure from the BG picture. In this approach the usual BG exponent, 
exp(— X/T), is deformed into the so called g-exponent (Tsallis distribution), 
exp g (— X/T) = [1 — (1 — q)X/T\ 1 ^ 1 ~ q \ such that for q — > 1 one recovers 
the BG picture again. 

The applications of the nonextensive statistical mechanics to nuclear 
and particle physics are numerous and we refer to [6] for details. In what 
follows we shall present the nonextensive version of the NJL model, the 
g-NJL model [7] q Details of the (/-NJL model and most of results were 
already presented in [7J, here we shall concentrate on the influence of dy- 
namical factors causing nonextensivity and represented by parameter q on 
the vicinity of the critical end point (CEP). 

2. Results 

Let us present first the basic elements of the g-NJL model introduced 
in [7] (to which we refer for more details). It is a q- version of standard 
SU(3) NJL model with U{1)a symmetry described in [I], with the usual 
Lagrangian of the NJL model used in a form suitable for the bosonization 
procedure (with four quarks interactions only), from which we obtain the 
gap equations for the constituent quark masses M^. 

Mi=mi- 2g s - 2g D {qjqj)(q k q k ) , (1) 

with cyclic permutation of k = u,d,s and with the quark condensates 
given by = —iTr[Si(p)] (Si(p) is the quark Green function); rrtj denotes 

the current mass of quark of flavor i. We consider a system of volume V, 
temperature T and the i th quark chemical potential \i\ characterized by the 
baryonic thermodynamic potential of the grand canonical ensemble (with 
quark density equal to pi = Ni/V , the baryonic chemical potential pb = 
\{p u + Pd + Ps) and the baryonic matter density as ps = \{p u + Pd + Ps)), 

fl(T,V,p l )=E-TS- ^ N i- ( 2 ) 

i=u,d,s 

2 Nonextensive calculations using the Walecka model for dense nuclear matter has 
been done in [8]. 
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The internal energy, E, the entropy, S, and the particle number, Ni, are 
given by |4j (here Ei - 
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The 77j = 1/ {exp [/3 (Ei - fa)] + 1} and fii = 1/ {exp [(@(Ei + /ij)] + 1} are, 
respectively, quark and antiquark occupation numbers with which one cal- 
culates values of the quark condensates present in Eq. ([I]) , 
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Eqs. ([I]) and (0) form a self consistent set of equations from which one 
gets the effective quark masses Mi and values of the corresponding quark 
condensates. 

The values of the pressure, P, and the energy density, e, are defined as: 



P(^,T) = -^^1, e( tM ,T) = ^^l with P(0,0) 



e(0,0) =0. (7) 



V ' ^" ' v 

The ^-statistics is introduced by using the q-form of quantum distributions 
for fermions (+1) and bosons (—1), namely, n q i = 1/ {e q ((3(Ei — fii)) ± 1]} 
where (x = 0(E - /x)) e q (x) = [1 + \(q - l)x\] x/l{q - 1)xl [7]. With such 
choice one can treat consistently on the same footing quarks and antiquarks, 
which should show the particle-hole symmetry observed in the g-Fermi dis- 
tribution in plasma containing both particles and antiparticles, namely that 
n q (E, (3, fx, q) = 1 — ri2- q (—E, (3, —fi) □. The g-NJL model is obtained by re- 



3 In a system containing both particles and antiparticles both q and 2 — q occur (i.e., 
one can encounter both q > 1 and q < 1 at the same time). It means that not only 
the q > 1 but also q < 1 (or (2 — q) > 1 have physical meaning in the systems we 
are considering what differs our q-NJL model from the q- version of the QHD-I model 
presented in [3- Notice that for q — > 1 one recovers the standard FD distribution, 
n(fi, T). It is important to realize that for T — > one always gets n q (fj,, T) — > n(u, T), 
irrespectively of the value of q [8], i.e., we can expect any nonextensive signature only 
for high enough temperatures. 



4 



MAZURIAN' JRGW PRINTED ON NOVEMBER 16, 2009 




Fig. 1. The pressure at critical temperature T cr as a function of compression p/po 
calculated for different values of the nonextensivity parameter q (the area marked 
at the left panel is shown in detail at the right panel). The dots indicate positions of 
the inflection points for which first derivative of pressure in compression vanishes. 
As in [4 for q = 1 the corresponding compression is p/po — 1.67 (and this leads to 
p = 318.5 MeV); it remains the same for q > 1 considered here (but now p = 321 
MeV for q = 1.01 and p = 326.1 MeV for q = 1.02) whereas it is shifted to 
p/p = 1.72 for q < 1 (p = 313 MeV for q = 0.99 and p = 307.7 MeV for q = 0.98). 



placing the formulas of Section [T] with their g-counterparts in what concerns 
the form of the FD distributions. Additionally, when calculating energies 
and condensates we follow [9] and use the g-versions of energies and quark 
condensates replacing Eqs. ([3]) and © by: 
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On the other hand, again following [9], densities which are given by the 
the q- version of Eq. © are calculated with n 9 's (not with n^, as in ([8]) and 
in ©). The pressure for given q is calculated using the above E q and the 
g-entropy version of Eq. Q with (cf . [TU] ) 



n q qi ln 9 n qi + (1 - n qi ) q ln 9 (l 



n, 



+ {n qi l-n q i} . (10) 



mazurian'JRGW printed on November 16, 2009 



5 




31 3 31 7 321 

M(MeV) 



2.4 
2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.! 



q=1 









17 



318 319 
u(MeV) 



320 



2.4 
2.2 
2.0 
1.8 
1.6 
1.4 
1.2 

H 

2.4 
2.2 
2.0 
1.8 
1.6 
1.4 
1.2 

H. 



q=0.98 


. . • • " 




** 







06 



307 308 
u(MeV) 



309 



q=1.02 


9 » • • ■ * 







24 



325 



326 
u(MeV) 



327 



328 



Fig. 2. The baryon compression p/po (calculated in the vicinity of the critical 
values of temperature and density indicated by the corresponding dotted lines) 
as function of the chemical potential p for different values of the nonextensivity 
parameter, q — 0.98, 1.00, 1.02. The summary presented in the top- left panel is 
detailed in the three consecutive panels. 



As example of how such approach works we present in Fig. [T]the pressure 
at critical temperature T cr as a function of compression p/po calculated for 
different values of the nonextensivity parameter q (see [7j for more details). 
We see that for q < 1 the critical pressure is smaller but for q > 1 it is bigger 
than the critical pressure for BG distribution. According to [7] it is directly 
connected to specific correlation for q < 1 and fluctuations for q > 1. The 
role of these factors is shown in more detail in Fig. [2 Notice the remarkable 
difference for the density derivative at the critical point: from the smooth 
transition through the critical point for q < 1 to a big jump in density for 
critical value of chemical potential for q > 1. It reflects the infinite values 
of the baryon number susceptibility, xb- 

xb = ^ E {dfH/d» B ) T = ~ E 9 2 n/d 2 p B \ T . (11) 

i=u,d,s i=u,d,s 

The transition between confined and deconfined phases and/or chiral phase 
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transition pQ can be seen by measuring, event by event, the difference in 
the magnitude of local fluctuation of the net baryon number in heavy ion 
collision They are initiated and driven mainly by the quark number 

fluctuation, described here by xb, and can survive through the freezout [TTJ 
I12j . Consequently, our q-NJL model allows to make the fine tuning for the 
magnitude of baryon number fluctuations (measured, for example, by the 
charge fluctuations of protons) and to find the characteristic for this system 
value of the parameter q. However, it does not allow to differentiate between 
possible dynamical mechanisms of baryon fluctuation. We close by noticing 
that using q dependent xb leads to g-dependent parameter e of the critical 
exponents which describe the behavior of baryon number susceptibilities 
near the critical point [13]. Whereas in the mean field universality class one 
has e = e' = 2/3 our preliminary results using q-NJL model show smaller 
value of this parameter for q > 1, (e ~ 0.6 for q=1.02) and greater for q < 1 
(e ~ 0.8 for q=0.98). It would be interesting to deduce the corresponding 
values of q from different models and from results on lattice with different 
universality classes - we plan to present it elsewhere. 
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